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We demonstrate that current and planned underground neutrino experiments could offer a pow¬ 
erful probe of few-MeV dark matter when combined with a nearby high-intensity low-to-medium 
energy electron accelerator. This experimental setup, an underground beam-dump experiment, is 
capable of decisively testing the thermal freeze-out mechanism for several natural dark matter sce¬ 
narios in this mass range. We present the sensitivity reach in terms of the mass-coupling parameter 
space of existing and planned detectors, such as Super-K, SNO-I-, and JUNO, in conjunction with 
a hypothetical 100 MeV energy accelerator. This setup can also greatly extend the sensitivity of 
direct searches for new light weakly-coupled force-carriers independently of their connection to dark 
matter. 


I. INTRODUCTION 


The existence of Dark Matter (DM) is clear evidence 
of physics beyond Standard Model (SM) and has inspired 
a major experimental effort to to uncover its particle na¬ 
ture. If DM achieves thermal equilibrium with the SM in 
the early universe, its present-day abundance can arise 
from DM annihilation with characteristic cross section 
CTU ~ 3 X 10“^®cm^/s. Alternatively if its abundance 
at late times is set by a primordial particle-antiparticle 
asymmetry, a thermal origin requires at least as large of 
an annihilation rate to avoid cosmological overproduc¬ 
tion. For either scenario, this requirement sets a pre¬ 
dictive target of opportunity to search for many of the 
simplest light DM models. 

Current and planned direct and indirect detection, 
and collider experiments will cover a vast subset of DM 
masses motivated by the thermal origin paradigm. How¬ 
ever, signihcant gaps remain in our current search strate¬ 
gies for low-mass DM. Indeed, the MeV-to-GeV DM mass 
window remains an elusive blind spot in the current 
search effort [T], despite the existence of viable mod¬ 
els EHS] - including those invoked to explain the ex¬ 
cess 511 keV photon line from the galactic bulge [9] 
with MeV scale DM annihilation into electron-positron 
pairs [1 g. Recent progress in our understanding of 
the status of MeV-scale DM has come from a combi¬ 
nation of re-interpretation of surface-level proton-beam 
neutrino experiments results HDHS], rare meson decays 
[HHni, electron beam dump experiments B- 

factories nana, precision measurements |51[ini[21], the 
CMB [251129) . and DM-electron scattering in direct de¬ 
tection experiments [30] . 

In this paper we propose a powerful new setup depicted 
schematically in Fig.[2 — the combination of a large un¬ 
derground detector such as those housed in neutrino fa¬ 
cilities and a low-energy (10-100 MeV) but high inten¬ 
sity electron-beam — which is capable of sharply testing 
the thermal origin scenario below ^ few 10s of MeV in 
DM mass. While our proposal requires a substantial in¬ 
vestment in an accelerator and beam-dump deep under¬ 
ground, it can signihcantly surpass the sensitivity of all 



FIG. 1. Schematic diagram of the proposed setup: a high 
intensity electron accelerator is placed in the vicinity of a 
large, underground neutrino detector. The electron beam im¬ 
pinges on a fixed-target or beam-dump to produce a dark 
force-carrier A', which can decay either visibly to e''"e“ or to 
DM particles. If the A' decays visibly and is long lived, it can 
enter the detector and directly deposit a large electromag¬ 
netic signal. If the A' decays invisibly to DM, the daughter 
particles inherit forward-peaked kinematics and scatter in the 
detector inducing observable target-particle recoils. 

other existing efforts in this mass range. This concept 
complements the DAE(5ALUS light-DM proposal |3T] in¬ 
volving an underground proton beam as well as other 
underground accelerator concepts [32H34] with different 
physics goals. 

For concreteness, we consider light DM that interacts 
with the visible sector through a kinetically-mixed |35j 
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FIG. 2. Sensitivity production for 10^^ electrons with 100-MeV energies impinging on an aluminum target positioned 10 m 
near the SNO+, JUNO, and SuperK detectors - since the latter two have comparable fiducial volumes, their projections are 
presented as a common curve. We conservatively assume thresholds of Er > 10 MeV for which the backgrounds are negligible. 
The CMB and direct detection constraints assume x/v’ constitutes all of the dark matter and regions above the relic curve 
correspond to parameter space for which each scenario can accommodate a subdominant fraction of the total DM (note that 
for subdominant DM, the CMB and direct detection bounds would also weaken). For the pseudo-Dirac scenario the relic 
curve was computed assuming only a small mass-splitting (100 keV < A < between the states that couple to the A! 

so standard freeze out is largely unaffected, but scattering at direct detection experiments is kinematically inaccessible. Since 
Kaon, mono-photon, and beam-dump constraints don’t scale as y, we conservatively adopt old = 0.5 and = 3 to 

not overstate these bounds; the plotted arrows show how the constraint moves when the product OLuirn^lmA'^^ is reduced by 
a factor of ten. The dotted LSND x SIDM curve denotes where the LSND bound shifts if old is chosen to satisfy the bound on 
DM self interactions ^0.1 cm^/g instead of the nominal aD = 0.5 which is conservative in other regions of parameter 

space. Note that for scalar inelastic DM, the key difference relative to the right panel is that the XenonlO region disappears as 
the scattering can be turned off. 


massive dark-photon A! [35]. Since light DM requires 
a comparably light mediator to avoid overclosure, this 
starting point loses no essential generality and our re¬ 
sults are qualitatively similar for different mediators. The 
most general renormalizable Lagrangian for this dark sec¬ 
tor contains 

•^D D H-+ Fdm, (1) 

where A' is the dark photon that mediates an abelian 
U{1)d force, and are 

the dark and hypercharge field strength tensors, and 
are the appropriate dark sector masses. After elec- 
troweak symmetry breaking, diagonalizing the gauge bo¬ 
son mass matrices induces a kinetic mixing with the pho¬ 
ton field strength e = ey cos 0\y , where 0w is the weak 
mixing angle. The DM Lagrangian contains a fermionic 
or bosonic MeV-scale DM particle charged under U{1)d, 

x{ip — w^)x, fermionic DM, , . 
\D^ip\'^ — bosonic DM, ^ ' 

where = 9^ + ig'A'^ is the covariant derivative. These 
simplest realizations of Cd assume a Dirac fermion or 
complex scalar DM states, but the model can readily 


be generalized to the “split” states of Majorana/pseudo- 
Dirac fermions or real scalars, in which case A! can cou¬ 
ple off-diagonally (inelastically) to the different mass- 
eigenstates and should be understood as a ma¬ 

trix acting on the split states. Moreover, each variation 
above can be particle/antiparticle asymmetric, which al¬ 
lows for weaker bounds from late-time annihilations into 
the CMB than the symmetric case [26j . 

One of the most important questions for such a model 
is the hierarchy of masses Ld- If rnA' < tTixlv^ 
mediator is the lightest state in the dark sector, so it will 
decay into SM particles. In this regime, the annihila¬ 
tion process that sets the relic density is t-channel (e.g. 
XX —f A'A') and, thus, independent of the mediator’s 
coupling to the SM. However, if niA' > ffis^i tfi® 

relic density is achieved through XX e^e~ annihila¬ 
tion, which proceeds via a virtual s-channel A' exchange 
and depends on both DM and SM couplings to the medi- 
atoiQ This latter scenario is predictive: since dark sec¬ 
tor couplings are bounded by perturbativity, sufficient 


^ In a certain region of parameter space, the sce¬ 

nario can still achieve the observed relic abundance through 
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FIG. 3. Same parametrization and conventions as Fig. but 
for asymmetric Dirac fermion DM; the symmetric version of 
this scenario is ruled out by the CMB. Note that the relic line 
has been replaced with the /otyer bound from the CMB. Since 
the residual antiparticle density in the early universe scales 
as ~ exp(—cru), a larger annihilation rate exponentially 
weakens the CMB constraint. Unlike Fig. all unshaded 
points are consistent with achieving the full DM abundance. 
For scalar asymmetric DM (not shown), the plot is similar, 
only the CMB bound is weaker and the “overclosure” curve 
is shifted upwards to match the “relic density” curve in Fig.[^ 
(right). 


experimentally sensitivity to the SM-mediator coupling 
can discoverer or decisively rule out this class of direct- 
annihilation models. 

The leading s-wave annihilation cross section for Dirac 
DM is given by 


av 


l6TTaaD'm,^{l + 

{m\, — )2 



(3) 


while for the scalar case the annihilation is p-wave 


av 


Stt e^aaDm'^v‘^{l -I- TOg/2m^) 
3 {m\, - 4m2 )2 


Y^l - ml/m; 


(4) 


and, hence, suppressed at late times (e.g. during CMB 
last scattering). Here an = and v is the rel¬ 

ative velocity between annihilating particles. These ex¬ 
pressions are approximately valid for the inelastic varia¬ 
tions of each scenario in the limit where the mass splitting 
between annihilation partners is small compared to their 
mass. Equating this annihilation rate to ^ 1 pbn x c 


XX A'A' annihilation in the “forbidden” channel [33, but 

this possibility is beyond the scope of our work. 


gives an important relation among the parameters of the 
model that ensures the correct cosmological abundance. 

This class of models can easily be adapted to cover a 
wide range of variations that can arise in a generic dark 
sector with a thermal history. Beyond the spin of the DM 
and mediator, we can classify the dominant dark species 
by its present day abundance (particle-antiparticle sym¬ 
metric vs. asymmetric) and the nature of its coupling to 
the mediator (elastic vs. inelastic). Since the experimen¬ 
tal setup we advocate in this paper (direct accelerator 
production and downstream scattering) is largely insen¬ 
sitive to these differences, we regard our approach as a 
convenient simplified model; the sensitivity of this ap¬ 
proach can be mapped onto most thermal DM scenarios 
in the low mass range. 

In comparing these scenarios against the thermal relic 
target, it is useful to adopt the ansatz introduced in [8] 
and define the dimensionless interaction strength 

/ \ 4 


_ 2 ( '^xlv ) 

y = e an \ - , 

V "lA' / 


(5) 


which, up to small additive corrections, is proportional 
to the annihilation rate and insensitive to assumptions 
about ratios of individual Largangian parameters. Thus, 
plotting y vs. invariantly presents the thermal 

reach of various experimental bounds. Some experi¬ 
ments, like direct detection, allow for unambiguous com¬ 
parison with the thermal target since the DM-SM non rel¬ 
ativistic scattering cross section is typically proportional 
to y. However, other constraints only depend only on a 
subset of the parameters that constitute y (e.g. collider 
production) so it’s necessary to make conservative as¬ 
sumptions about the other parameters for a robust com¬ 
parison with the thermal target - for a discussion see 
Sec. El 

For any DM search involving accelerator production 
and downstream detection, the signal yield is propor¬ 
tional to the number of DM particles produced and the 
combined probability that the DM reaches the detector 
and scatters off a stationary target inside. If the media¬ 
tor {A') can kinematically decay to DM then the signal 
yield scales as 


Signal CX POT X (PQeom.Acc.) X (-^Atonis in det. ),(6) 

where POT stands for the number of particles-on-target, 
P Geom. Acc. IS the probability that the produced DM path 
crosses the geometry of the detector, and iVAtoms in det. 
is the number of target particles inside such a detec¬ 
tor. The ideal light DM beam-dump experiment seeks 
to maximize the product of these factors, while minimiz¬ 
ing background rates. 

Our proposal, which pairs a high-intensity low-energy 
electron-beam accelerator with a large underground de¬ 
tector, Fig. E has several key features that maximize 
signal and reduce potential backgrounds: 


1. The next generation of deep underground neutrino 
detectors are among the largest detectors ever built 
or proposed. 
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2. Underground facilities, typically situated inside a 
mine at a few water-equivalent km underground, 
boast enormous overburdens, which exponentially 
suppress cosmic and other environmental back¬ 
grounds. 

3. An electron beam can be advantageous over a pro¬ 
ton beam for probing MeV-scale New Physics (NP) 
since beam-related backgrounds are significantly 
suppressed compared to the latter, as the produc¬ 
tion of neutrinos and neutrons is significantly less 
in the electron beam case. 

4. For DM masses small compared to the beam en¬ 
ergy, radiative production in electron-nucleus col¬ 
lisions is extremely forward-peaked and offers ex¬ 
cellent (order-one) geometric acceptance for fixed- 
targets placed in the vicinity of large neutrino de¬ 
tectors. 

5. A high-intensity electron-beam — e.g. a continuous 
wave (CW) beam at existing facilities — capable 
of running at mA currents combined with a large 
neutrino detector offers a significant improvement 
in luminosity compared to previous experiments. 

6. Finally, the proposed light DM search scheme does 
not interfere with other physics goals of large un¬ 
derground detectors. Moreover, for some facilities, 
a project of this scope and ambition may offer an 
entirely novel physics program - especially for de¬ 
tectors whose main physics goals have already been 
met (as is the case for most solar neutrino detec¬ 
tors). 


For the remainder of this paper we narrow our scope 
to simple models in Eq. Q, in which DM communicates 
with the SM through the vector portal. In Sec. [ll] we 
discuss the production and detection of a lightly cou¬ 
pled new particle at underground facilities; Sec. m dis¬ 
cusses the potential promising facilities for realizing the 
proposal in this paper; in Sec. IV we discuss the main 
backgrounds for such an experimental setup; In Sec. 
we discuss existing experimental constraints on this pa¬ 
rameter space and the projected sensitivity of our setup 
for various detectors; finally in Sec. VI we offer some 
concluding remarks. 


II. PRODUCTION AND DETECTION 
A. Invisibly Decaying Mediator 

In this section we follow a well-established routine for 
calculating the elastic scattering of light DM y produced 
from the decay of on-shell mediators A' [191 [20]. The pro¬ 
duction of mediators occurs via radiative emission from 
electron-nucleus scattering, Z + e~ Z + e~ + A', and 


Visibly Decaying A' 



FIG. 4. Constraints and projections for the visibly decaying 
A' -A scenario. The shaded region represent existing 

bounds from beam dump and accelerator searches taken from 
m- The red curve is the 10 event yield projection for the 
same JUNO or SnperK setnp described above, but deriving 
a signal yield from both visible decays and A' absorption. 
Comparable projections are obtained for other detectors con¬ 
sidered above. The dashed blue curve represents projections 
for future visibly-decaying A' searches reviewed in |4T] and 
the dashed black region represents model-dependent interpre¬ 
tations of supernova observations. If additional interactions 
allow the A' to thermalize with (heavier) DM states also pro¬ 
duced in during the explosion, the energy loss bounds can be 
alleviated- see discussion in Sec. m 


is calculable with a minimum of uncertainty, by stan¬ 
dard QED methods. For the results in Sec. |V| we cal¬ 
culate the production of DM at a fixed target using a 
modified version of Madgraph from Ref. [331 j which 
includes electron-nucleus scattering with form-factors ob¬ 
tained from Ref. m- 

The prompt decay of A! mediators to light DM, A -A 
XX, creates a flux of (semi-)relativistic x particles that 
subtends the detector volume. For invisibly decaying me¬ 
diators, the signal consists of detector particle recoils in¬ 
duced by A daughter particles - e.g. ye ye mediated 
by virtual A exchange. The scattering yield for this pro¬ 
cess is 




NA'Tiee 

Ns 



das 

dEn 


e{eD-e ,),(?) 


where Eq is the e~ beam energy, Ec is the recoil “cut” 
threshold, Na' is the total number of A produced, Ei is 
the daughter DM energy, E^ is the target’s recoil energy, 
Uf. is the target (electron) number density, and is the 
DM path length through the detector. For our signals of 
interest, typically E^ > 10 MeV, the detection efficiency 
e is of order-one. The theta function enforces geometric 
acceptance by admitting only particles whose angle 9^ 
with respect to the beam line falls inside do, the angle 
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subtended by the detector. We have presented this rate 
as an average over Ng simulated monte carlo events of A! 
production and decay in the target. The nature of the 
scattering cross section can vary depending on energy and 
detection methods. For most large neutrino detectors, 
the dominant signal is from electron recoils. Assuming 
Dirac fermion DM, the recoil spectrum is 

dos _ 87rQ^me(l - E^jEj) 

dER a ^ {m\, + 2meEji)'^ 

in the Ei ^ limit; a parametrically similar expression 
can be obtained for inelastic or scalar DM variations [20] . 
Other signatures may include elastic and quasi-elastic 
scattering of x on nucleons, and nuclear elastic recoils. 
For the inelastic-DM scenario, the process yie —>■ X 2 S- 
is well-approximated by the above formula in the limit 
where the mass splitting is small compared to the inci¬ 
dent DM energy. 


B. Visibly Decaying Mediator 


For the visibly decaying mediator (m^' < the sig¬ 

nal consists of both decays and absorptions. For a sim¬ 
ulated population of Na' mediator particles, the total 
yield is V = Ya + To, where the effective absorption {A) 
and decay (D) contributions are 


N. 


Ya = 


Yd = 


Ng 


e-^ ''dE^{E,)eMdD - 0,) 


NA'UeE 


2=1 


/E^ 


N. 


Ns 

i£ t-A rd; / rf; s 

— j 0(6 »d - fli) . 


(9) 


2=1 


Unlike in Eq. Q, here Ei/9i are now the A' energy/angle 
of the ***' A' produced in the simulation, Xi = (1 — 
is the A' boost factor, di is the path traveled between 
target and detector, and the visible A! width is 

r = ^ (1 + - Amllm\, . (10) 


The absorption cross section for the A'e —>■ ye process 
with a stationary detector electron is 


daA _ m^{\A\^) 
dEn 327r|p*ps ’ 


( 11 ) 


where Er is the electron recoil energy, s = {me + Ei)^, 
= [(s — TOg — — 4TOgm^,]/4s is the CM three- 

momentum in the lab frame, and the absorption ampli¬ 
tude A = Ag +At is a sum of s and t channel amplitudes 

^t= , ^ ^M Pnh^il^T -^7 + rneh''u{pT)£^'el, ( 12 ) 

(t me) 

■^^=r~Z-^'d{pR)x^{j/T + I^A' + meh''u{pT)£f el, (13) 


where Pt/r are the four-momenta of the target/recoil e , 
Pa' is the A' four-momentum, and e^ are gauge boson 
polarization vectors. 


III. SUITABLE FACILITIES 

The experimental strategy that we propose relies on 
production of light DM by an electron beam impinging on 
a target or beam dump. Further detection downstream 
is realized at a large detector. The advantage of this ap¬ 
proach is that it permits use of an existing or planned 
detector. We now discuss existing and upcoming detect¬ 
ing facilities which could be suitable to meet the physics 
targets outlined in this paper. 


A. Large water-based Cerenkov detectors 
Super-Kamiokande 

The Super-Kamiokande (SuperK) detector |42| in the 
Kamioka mine in Japan contains a 50,000-ton water 
Cerenkov detector, at a depth of 2700 m.w.e.. It was 
designed to measure neutrino rates from different sources 
such as those from astrophysics and the atmosphere; and 
for proton decay. It consists of a welded stainless-steel 
tank 39 m in diameter, it stands 42 m tall, and is sen¬ 
sitive to electron recoils above 4.5 MeV. An important 
feature of this type of detector is directional sensitivity, 
which can be used to further suppress environmental 
backgrounds. Super-K is the largest among the existing 
detectors suitable for our proposal. It currently hosts a 
very low-power medical-type electron linac for calibra¬ 
tion purposes [15] . 


Hyper-Kamiokande 

The next large water Cerenkov detector at Kamioka mine 
is the Hyper Kamiokande project [44], intended to also 
measure neutrino rates, look for proton decay, in addi¬ 
tion to being the far detector for a long baseline neu¬ 
trino oscillation experiment using neutrinos originated 
from the upgraded J-PARC. Its design is based on the 
same detection technology that SuperK featured. It has 
a 40m X 54m x 250m detector volume, with a hducial 
volume ^ 20 times greater than that of SuperK, located 
at a depth of 1750 m.w.e.. 


B. Large liquid scintillator detectors 
Borexino 

The borexino detector is located in the Gran Sasso 
laboratory in Italy, under 3800 m.w.e. and contains a 
large volume liquid scintillator in a 8.5 m diameter vessel 
[45] . The detector registers energy depositions (mainly 
y-events and electron recoils, and with some quenching 
factor penalty it can also observe proton recoils) and 
offers good energy resolution with a low threshold (^200 
keV), which allowed Borexino to measure different 
components in the solar neutrino flux. It is cur¬ 
rently being considered as a short-baseline detector for 
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neutrinos emitted from an intense radioactive source |46| . 

Kamland 

The Kamland detector is located in Kamioka mine, a 
13 m sphere filled with liquid scintillator. It has played 
an important role in the precision determination of 
neutrino mass splitting and mixing [47] , It currently has 
been re-profiled as a tool to search for the neutrino-less 
double-beta decay [iS] , 

SNO+ 

SNO-I-, the upgraded Sudbury Neutrino Observatory 
(SNO) is a 12 m diameter liquid scintillator, situated at 
« 6000 m.w.e. [13]. The detector features a significant 
light-yield upgrade over its predecessor, and chief 
amongst its current science goals are neutrino-flux mea¬ 
surements from the pp, CNO cycle, from supernovae, 
in addition to searching for neutrino-less double-beta 
decay. SNO-I- energy thresholds are expected to be near 
200 keV. 

JUNO 

The Jiangmen Underground Neutrino Observatory 
(JUNO) is a multi-purpose neutrino oscillation exper¬ 
iment, m situated 700 m underground and currently 
under construction, with data-taking expected to com¬ 
mence in 2020. The detector itself is a 35.5 m diameter 
liquid scintillator embedded in a water pool which will 
also serve as a Cerenkov muon veto. 

DUNE 

Another interesting opportunity may arise from the 
DUNE(LBNF) project. A detector containing 40 kton 
of liquid argon is planned to be built at the Sanford Un¬ 
derground Research Facility (SURF), at a depth of 4300 
m w.e. While the main goal of this detector will be ob¬ 
serving long baseline neutrino oscillations, it can also be 
used in conjunction with a linac at SURF. 


C. Large-scale dark matter detectors and Qv2l3 
detectors 

All of the detectors described above have over 100 tons 
(sometimes several kilotons) of detecting material. It is 
also worth noting that underground facilities often house 
very large, 0(1 ton), detectors aimed at the direct dark 
matter detection and/or neutrino-less double-beta de¬ 
cay processes. These include XFNONIT [^, LUX [ 52 ] . 
DEAP-3600 [ 55 ] , future LZ experiment [ 53 ] - all devoted 
to DM direct-detection for > 10 GeV masses - as well 
as other experiments such as e.g. EXO [ 55 ] and CUORE 
[ 55 ] that take aim at 0v2j3. Although the active mass in 
such detectors is smaller than neutrino detectors, some 
of them have remarkably low energy thresholds, and are 
capable of detecting coherent scattering on nuclei. This 
coherence typically adds another enhancement factor of 


Z ~ 0(50), compared to incoherent scattering on elec¬ 
trons, which helps to offset their small mass. In many of 
the underground facilities these detectors co-exist with 
large neutrino detectors, and therefore it may be possi¬ 
ble to perform an accelerator-based search for light DM 
using multiple detectors hosted in the same laboratory. 


IV. BACKGROUNDS FOR UNDERGROUND 
ACGELERATOR BEAM DUMP SEARGHES 

One of the advantages of the experimental setup pro¬ 
posed in Sec. |ll|is the clean environment that the under¬ 
ground facilities offer. In this section we discuss the rele¬ 
vant backgrounds that give rise to SM muons and neutri¬ 
nos, which can mimic the signal from light DM scattering 
off electrons. We split the discussion into beam-related 
backgrounds arising from reactions in the target or dump 
and environmental backgrounds arising from ambient ra¬ 
dioactivity and cosmic-rays traversing the detector. 

We first discuss beam-related backgrounds. For the 
proposal in this study which relies on an electron beam 
running at an energy below beam-related back¬ 

grounds may arise from a few sources. Direct production 
of neutrinos in reactions due to weak neutral or charged 
currents, eN eNvv and eN —>■ v^^N' + X, followed 
by neutrino scattering in a detector, vN —?> vN' -|- A or 
ve —>■ ve, is small. Indeed, the production cross section is 
« /6/nucleon, [57], giving only up to 10® — 10® neutrinos 
with El, ~ beam energy, produced in one year’s running 
and assuming an aluminum beam-dump. For a neutrino 
to fake the DM signal, it needs to subsequently scatter in 
the detector situated behind the dump. The probability 
for a neutrino to scatter off an electron in the detector 
is given by where no and Id are the num¬ 

ber density and longitudinal length of the detector. The 
scattering cross section scales as GpE„me. Even if the 
neutrino produced in the target carries most of the beam 
energy, the event rate from this type of background will 
be well below the atmospheric neutrino background. 

A more serious source of beam-related backgrounds is 
the spallation of nucleons by the electron beam. This 
may lead to a neutron scattering background and to Gp- 
unsuppressed production of neutrinos by resulting 
emitters. The inelastic cross sections of electrons on 
nuclei can reach ^ mb range, and therefore, the num¬ 
ber of produced neutrons and neutrinos will significantly 
exceed the direct production discussed in the previous 
paragraph. However, extra neutrons produced this way 
will be intercepted by some meters of rock separating 
the production point from the detector. The neutrino 
energies will be limited to Ui, < 5 MeV, typical for /3- 
processes, and can be removed by a simple recoil energy 
cut. Therefore, both types of backgrounds are reducible. 
Moreover, a judicial choice of beam dump material can 
reduce this contamination even further, so we conclude 
that the beam-related backgrounds for the experimental 
setup we advocate for are negligible. 
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The nature of the accelerator beam plays a crucial role 
in whether environmental backgrounds are important or 
not. The higest intensity beams may only be attained by 
using a continuous wave (CW) beam. In that scenario, 
given the 0{1) duty cycles common to a CW beam, us¬ 
ing timing to reduce cosmic-related backgrounds is not 
an available handle, unlike at a pulsed beam. Never¬ 
theless, environmental backgrounds can still be reduced 
to a negligible level through a combination of existing 
techniques. We base the following feasibility argument 
on the lessons learned from the results by the Borex- 
ino collaboration. First, the radiated massive A' from 
the reaction eN —>■ eNA' will carry most of the beam 
energy. Assuming a 100 MeV electron beam, this means 
that one could use a high electron recoil-energy threshold 
in the detector. In particular, a En > 5MeV threshold 
would eliminate all radiogenic backgrounds that Borex- 
ino has encountered thus far. This still leaves cosmic ray 
muons, cosmogenic neutrons, and solar neutrinos origi¬ 
nating from the + Ve reaction in the sun as a 

potential source of contamination for a light DM scatter¬ 
ing signal. However, these could be significantly reduced 
through a combination of shielding, fiducialization, tim¬ 
ing veto following activity in the detector (for stopped 
muons), and directionality, as was attained in Ref. [58] 
(see Fig. 5). In Sec.j^we conservatively adopt a 10 MeV 
energy threshold, and assume a background-free experi¬ 
mental setup for our sensitivity estimates. 


V. PROJECTED SENSITIVITY AND RESULTS 

We now turn our attention to the existing constraints 
on DM in the (1-100) MeV mass range for the class of 
models in Eq. 0- One of the strongest limits for DM 
masses below ttIt^ comes from a reinterpretation of re¬ 
sults from surface-level neutrino experiments, in partic¬ 
ular, LSND’s measurement of neutrino-electron neutral 
current scattering and recast by m- Production of 
DM in such a setup proceeds from protons impinging on 
a fixed target. The latter gives rise to mesons in the fi¬ 
nal state, primarily tt^s which consequently can decay to 
either a DM pair or to a mediator via 7r° —^ 
or 7r° —)■ yA', and finally with decay of the A' into a 
DM pair. Once produced, DM detection proceeds anal¬ 
ogously to neutrinos, with DM traveling through a near 
or far detector that is sensitive to the scattering process 
xT —)■ yT. Here T could be a nucleus, nucleon, or an 
electron. The results above use the electron-scattering 
channel. A proposal for further running at MiniBoone 
could further improve the constraints at neutrino exper¬ 
iments for higher mass DM [60] . 

In the context of DM searches, one of the main lim¬ 
itations of neutrino factories is that neutrinos can be a 
background to the DM signal. In contrast, electron beams 
do not suffer from this limitation [19]. In fact, for an elec¬ 
tron beam with energy of 100 MeV, neutrinos originating 
from muons or pions are no longer kinematically accessi¬ 


ble. One of the strongest constraints on light DM comes 
in fact from E137 in a manner analogous to neutrino ex¬ 
periments m where DM production in the target pro¬ 
ceeds via eN —>• eN A' fxxlv’^* ■ E137 sets a constraint 
on DM from the non-observation of anomalous electron- 
recoil events in a detector ~ 300 meters downstream of 
the target. 

Another relevant constraint to the parameter space 
that we focus on is that from B-factories. In particular, 
Babar’s search for untagged T(3S') decays to x+invisible 
EH was reinterpreted by Refs. mills] to place limits on 
DM where the DM could be produced in the reaction 
e“''e“ —>■ yA' (and similarly for off-shell production of 
DM). Kaon decays, and in particular the missing energy 
mode K'^ —>■ also impose some consitraints on the 

parameter space of the model via ^ -k'^xxIh^E* ■ 

While standard nuclear direct detection experiments 
are not sensitive to DM masses near the MeV scale, elec¬ 
tron scattering does constrain masses this light. While 
these searches are background-limited, they still set a 
conservative limit. In particular. Ref. [30] reinterprets 
XENONIO electron data to set strong limits. 

Other constraints on the kinetic mixing scenario can 
be derived from LEP which accurately measured the Z 
boson mass. If the kinetic mixing parameter is too large 
it can induce an unobserved shift in mz^ which constrains 
e < 10-2 [62]. 

In Figs. and we present the expected sensitivity 
of the combination of an 100 MeV electron beam im¬ 
pinging on a beam-dump situated « 10 meters from a 
large underground detector, such as that of the SNO-f, 
JUNO, and Super-K experiments. These results are 
plotted as a function of the DM mass and the variable 
y = a£)e^(m^/mA')^ which, up to small additive correc¬ 
tions, is proportional to the annihilation rate. With this 
anzatz, the thermal relic target is invariant on the y vs. 

plane and we don’t need to assume any particular 
values of individual inputs for a robust comparison with 
experiment. 

However, some constraints on e (e.g. from LEP or 
kaon decays) don’t depend on au or on the mass ra¬ 
tio we must make conservative assumptions for 

these inputs in order to construct y for these bounds. 
Eor the predictive regime of interest, m^' > 
the thermal relic annihilation rate is proportional to 

ajorr?^!m\, ^ so we choose both ajo and the ratio 
my^/mA' to be 0(1), which reveals the maximum param¬ 
eter space on the y vs. my- plane; it is always trivially 
easy to overstate these bounds by choosing smaller values 
of these quantities, which shifts the y contours for col¬ 
lider and kaon bounds downwards by a large amount, so 
this approach is the most conservative. Note that Bullet- 
cluster constraints on DM self-interactions [63] [64] set an 
upper bound on ajo and therefore on y, which we include 
in Fig. [^ 

Figure [^ considers two scenarios: scalar thermal relic 
DM (right) and Majorana-like (or pseudo-Dirac) DM 


(left) with a mass-splitting A > 100 keV. The latter 
scenario is not subject to constraints from direct detec¬ 
tion - because the upscatter is kinematically forbidden at 
tree level - or from the CMB because the annihilation is 
exponentially suppressed when the heavier, unstable an¬ 
nihilation partner is depopulated before last-scattering. 
For scalar thermal relic DM, the minimal symmetric real¬ 
ization is cosmologically viable because the p-wave anni¬ 
hilation rate is sharply redshifted at the epoch of recom¬ 
bination, thereby greatly alleviating the CMB bound. 

In Fig. [^we show a similar plot for asymmetric Dirac 
fermion DM coupled to a dark photon. Unlike the 
particle-symmetric Dirac scenario (not shown), which is 
ruled out by the CMB for DM masses below 10 GeV 
[m [23111 [05], this variant remains viable because anni¬ 
hilations at early times exponentially deplete the antipar¬ 
ticle density before last scattering so the energy injected 
into the CMB at late times is reduced. Since the antipar¬ 
ticle density scales as ~ exp(—tru) oc exp(—p), larger 
values of y are less constrained, so the CMB imposes a 
lower bound on y for this class of models. In setting this 
limit, we demand that av/m^ ~ 10“^®cm^/s/GeV [H] 
while fixing the observed DM abundance at late times 

m- 

The projected reach of the proposed setup is quite sig¬ 
nificant. For DM lighter than 10s of MeV, our proposed 
production and re-scattering search can exceed the LSND 
reach by many orders of magnitude and completely cover 
all viable thermal scenarios whose freeze-out abundance 
is set by e^e“ annihilation. Only a dedi¬ 

cated missing momentum fixed-target scheme, Ref. [H], 
can offer more sensitivity since the yield for that setup 
scales as instead of e^. 

Finally, for completeness we also include the con¬ 
straints on the visibly-decaying A! scenario in Fig. 

We observe that the proposed setup would be the most 


sensitive direct probe of this parameter space, greatly 
improving the reach of previous beam dump constraints. 
Although this region is naively excluded by theorists’ in¬ 
terpretations of the supernova cooling bounds from SN- 
1987a, this constraint is model dependent, so we present 
it as a dashed black line in Fig. [^ For example, if in ad¬ 
dition to A there is an additional mediator A" and light 
DM X that are fully thermalized during the SN explo¬ 
sion, the A-x interaction can also keep A thermalized 
and quench the energy loss through A emission. There¬ 
fore, the sensitivity reach to visible decays adds value to 
the proposed set-up. 

VI. DISCUSSION 

This paper proposes a new experimental set-up to 
probe MeV-scale DM. The powerful combination of a 
low-energy high-intensity electron beam situated near 
a large scale detector underground can be complemen¬ 
tary to existing proposals using electron direct detection 
[301 [67] , fixed-target missing momentum [H] and future 
B-factories [Hj (for similar proposals using proton beams 
underground see Ref [32] and Ref. m)- In particular, 
it is capable of surpassing the sensitivity of LSND, a 
surface-level neutrino detector benefiting from a high- 
intensity proton beam, and decisively probing thermal 
DM lighter than ~ 10 s of MeV in some of the simplest 
light DM scenarios. Given that many underground de¬ 
tectors have already achieved their original physics goals, 
a light-DM search via accelerator production may add an 
exciting new application to these facilities. 
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